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Consider  as  e  *  O'* ,  the  solution 

y(x)  -  {1  -  exp( -x/e )}  /  {1  -  exp(-1/e)> 

of  the  following  singularly-perturbed  linear  two-point  boundary-value  problem 

e  y"(x)  +  y'(x)  *  0,  y( 0 )  -  0,  y(1)  -  1  . 

When  this  problem  is  cast  as  a  system  of  first-order  differential  equations, 
using  the  substitution  y^  =  y  and  y2  »  y' ,  one  is  led  to  consider  the 
problem 


i  n  a  fy«<*>l  n  .  y,(x) 

r 1  0i  a_  1  r0  1 i  1 
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o,  yi(1)! 


y2cu 


As  e  0+  the  solution  of  this  first  order  system  has  the  property  that 

y2(0)  blows  up  like  e”1  while  y^(x)  remains  uniformly  bounded  on  [0,1]. 

This  paper  presents  a  simple  constructive  method  of  solving  such  singularly- 
perturbed  linear  two-point  boundary-value  problems;  a  method  which  has  been 
generalized  to  analyze  the  solution  of  such  problems  by  finite  different 
schemes  based  on  Euler's  method. 
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SIGNIFICANCE  AND  EXPLANATION 


This  paper  presents  a  simple  constructive  proof  of  the  existence  of 
solutions  to  a  class  of  singularly-perturbed  linear  two-point  boundary-value 
problems.  Such  problems  arise  when  detailed  models  of  physical  phenomena 
involve  effects  that  occur  on  markedly  different  temporal  or  spatial  scales, 
the  (singularly)  small  parameter  measuring  the  disparity  between  the  scales. 
This  simple  constructive  proof,  through  slight  modification,  also  allows  one 
to  prove  the  existence  of  solutions  of  a  class  of  difference  schemes  used  to 
approximate  the  solution  of  the  continuous  problem. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


ANALYSIS  OF  A  SINGULARLY-PERTURBED  LINEAR 
TWO-POINT  BOUNDARY-VALUE  PROBLEM 


Warren  E.  Ferguson,  Jr* 

1*  Introduction 

When  a  high-order  singularly-perturbed  linear  two-point  boundary-value  problem  is  cast 
as  a  first-order  system  of  differential  equations,  it  is  usually  the  case  that  the  solution 
of  this  first  order  system  does  not  remain  uniformly  bounded  as  the  small  parameter  tends 
to  zero.  This  leads  us  to  consider,  as  e  +  0+,  a  differential  equation  on  [0,1]  of  the 
form 
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subject  to  appropriate  boundary  conditions.  Here  it  is  assumed  that  the  real  parts  of  the 
eigenvalues  of  A^(x)  are  negative  while  those  of  A33(x)  are  positive.  In  this  paper 
reasonable  general  conditions  are  placed  on  the  boundary  condition  which  ensure  that  the 
solution  y  of  this  differential  equation  satisfies  the  followlngt 


Requirement  1.1;  For  all  sufficiently  small  positive  e 


exists  and  is  finite. 


lim+ 

e+0 


ey^O)' 

y2(0) 

ey3(1) 


Problems  of  this  type  have  been  considered  by  other  authors  [6,11].  What 
distinguishes  the  analysis  presented  in  this  paper  is  the  fact  that  a  straightforward 
analog  can  also  be  used  to  analyze  the  stability  and  consistency  properties  of  difference 
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schemes  used  to  solve  this  problem. 

In  section  2  we  precisely  define  the  class  of  singularly-perturbed  linear  two-point 
boundary-value  problems  to  be  studied.  In  sections  3  and  4  a  simple  constructive  technique 
is  presented  which  allows  one  to  characterize,  in  a  reasonable  general  manner,  when  such 
singularly-perturbed  boundary-value  problems  admit  a  solution  satisfying  Requirement  1.1. 

In  section  5  applications  of  the  results  presented  in  sections  3  and  4  are  described. 
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2.  Preliminaries 

Consider  the  following  two  boundary-value  problems  BV(v)  for  v  =  G,  S  which  differ 
only  in  the  boundary  conditions  imposed  on  their  solutions : 

Ly  -  f  and  8^v'y  -  g*v*  t 

where 
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and 

8<V)v  ^  B<u)y(0>  +  8^v)y(1>  . 

The  unigue  solvability  of  the  general  problem  BV(G)  will  be  related  to  the  unigue 
solvability  of  the  special  problem  BV(S)  under  the  following  assumptions. 


assumptions  2.1:  For  all  sufficiently  small  positive  e : 


0) 

n  * 

n1  +  n2  +  n3 * 

1) 

the 

components  y^(x), 

fi(x),  and  g!v) 

of  the  n-vectors 

y(x) 

,  f(x),  and  g*V) 

respectively  are 

n^-vectors , 

2) 

the 

order  by 

matrices  A^j(x) 

»  A| jtxie )  depend 

smoothly  on  x  and  e 

/ 

3) 

the 

vector  f(x)  «  f(x»e)  depends  smoothly  on  x  and  e, 

4) 

the 

order  n  matrices 

*SG) -<”«•>. 

b<s>  -  b<g>e 

and 

order  n  vectors 

i-  -i 

g(v)  =  g(v)(c). 

depend  smoothly  on 

5)  B, 


(S) 


ei 

(S) 

"o 

I 

0_ 

and  B 1  * 

0 

ei 

are  order  n  matrices. 


6)  for  some  positive  constant  u  independent  of  e  the  eigenvalues  of 
A^(x)  have  real  part  less  than  -u  while  the  eigenvalues  of  A22(x) 
have  real  part  greater  than  +u . 

By  smooth  dependence  on  x  and/or  e  we  mean  that  all  possible  derivatives  with  respect 
to  x  and/or  e  exist  and  are  continuous. 
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Whan  one  formally  sets  e  »  0  the  differential  equation  Ly  -  f  reduces  to  the 
following  familiar  reduced  system  of  algebraic-differential  equations 

y^x)  -  -A^(x)[A12(x)y2(x)  +  f^x)]  , 

d  *  * 

y2<x)  “  A22(x)y2(x)  +  f2(x)  ,  and 

y3<x)  -  -a"3(x) [A32(x)y2(x)  +  f3<x)] 

where 

A22(x)  H  A22(x)  -  A21(x)A^(x)A12(x)  -  A23(x)A33(x)A32(x)  ,  and 

f*(x)  =  f2<x)  -  A21(x)A;](x)f1(x)  -  A23(x)A^(x)f3(x)  . 

When  assumptions  2. 1  hold,  standard  results  in  singular  perturbation  theory  suggest  that 
as  e  ♦  0+  the  solution  of  BV(S)  approximately  satisfies  this  reduced  system  of 
algebraic-differential  equations  in  each  fixed  open  subinterval  of  [0,1]. 

The  proof  of  the  unique  solvability  of  BV(S)  begins  by  casting  the  BV(S)  as  an 
equivalent  operator  equation  acting  on  an  appropriate  Banach  space.  To  cast  BV(S)  into 
an  equivalent  operator  equation  we  will  use  the  fundamental  solution  matrices  Y^(x,s), 
Y2(x,s) ,  and  Y3(x,s)  defined  as  follows: 

t  I^Y^x.s)  -  A11(x)Y1(x,S)  ,  Y1(S,S)  «  I  , 

|^Y2(x,s)  -  A22<x)Y2(x,s)  ,  Y2(s,s)  =*  I  ,  and 

6  feY3<X'8)  “  *33(x)Vx's)  ,  Y3(.,s>  -  I  . 

Assumptions  2.1  and  well-known  results  for  initial-value  problems  guarantee  the  existence 
of  these  fundamental  solution  matrices. 

The  results  presented  in  this  paper  depend  primarily  on  three  fundamental  theorems. 
The  first  theorem  relates  the  unique  solvability  of  BV(G)  to  the  unique  solvability  of 
BV(S). 
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Theorem  2.2;  Let  8V(S)  be  uniquely  solvable.  Define  Y^s^  to  he  the  fundamental 
solution  matrix  of  BV(S),  that  is 

i.Y<S)  -  o  and  8(s>Y!s>  -  I  . 

Then  BV(G)  is  uniquely  solvable  if  and  only  if  the  matrix  g[G)y(S)  is  nonsingular. 

This  theorem,  whose  proof  can  be  found  in  [8,9],  was  used  by  Keller  to  show  that  a 
difference  scheme  is  stable  and  consistent  for  a  uniquely  solvable  boundary-value  problem 
if  and  only  if  it  is  stable  and  consistent  for  the  related  initial-value  problem.  The 
proof  of  this  theorem  for  difference  equations  rest  on  Kron's  method  of  tearing  [12]  i  the 
boundary  conditions  of  the  boundary-value  problem  are  torn  from  the  matrix  describing  the 
difference  scheme  and  replaced  by  the  boundary  conditions  of  the  related  initial-value 
problem.  This  idea  of  replacing  one  set  of  boundary  conditions  by  another  underlies  the 
proof  that  BV(G)  is  uniouely  solvable. 

The  second  theorem,  whose  proof  can  be  found  in  [4] ,  orovides  needed  estimates  on  the 
size  of  the  fundamental  solution  matrices  Y y  and  Y?. 

Theorem  2 . 3 i  There  exist  positive  constants  K  and  X,  independent  of  e,  such  that 


Y, (x,s) 1 

< 

„  -X  ( x-s )  /e 

1C  e 

for 

0  < 

a  < 

X  <  1 

1 

m 

* 

m 

m 

Y3(x,a) t 

< 

K  a-X  (s-x)/e 

for 

0  < 

m 

X  < 

m 

8  <  1 

for  all  sufficiently  small  positive  e . 

tn  this  theorem  the  norm  | • |  used  is  any  fixed  n-vector  norm,  say  the  infinity 
vector  norm.  Note  that  the  condition  that  e  be  sufficiently  small  is  important  [1,2] t 
for  unless  e  is  sufficiently  small  |Y.j(x,a)!  need  not  decay  with  increasing  (x-s) 
even  though  the  eigenvalues  of  1 ( x)  have  uniformly  neqative  real  part. 

The  third  theorem,  whose  proof  can  be  found  in  [7],  states  that  operators  in  a  Banach 
space  sufficiently  near  the  identity  are  invertible.  The  Banach  space  Cn  used  in  this 
theorem  consists  of  the  space  of  n-vector  valued  functions  of  x,  continuous  on  [0,1], 
endowed  with  the  norm  l«l  defined  bys 

I f ( x ) I  5  max  I f ( x ) I 

[0,1] 
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Theorem  2.4:  Let 
IBI  <  1.  Then  I 

This  theorem  will 
both  BV(S)  and 


B  be  an  operator  on  the  Banach  space  Cn  which  admits  the  bound 

•B  is  invertible  with 

00  N  H+ 1 

(l-B)"1  =  l  B1,  and  l(I-B)"1  -  l  B%  <  • 

i“0  1=0  1 

be  used  to  obtain  a  generalized  power  series  expansion  of  the  solution  of 
BV(G) . 


3.  Solution  of  the  Special  Problem 

Let  us  now  show  that,  for  all  sufficiently  small  positive  e ,  the  speical  problem 
9V(S)  is  uniquely  solvable.  ka  stated  in  section  2,  the  first  step  is  to  cast  BV(S)  as 
an  equivalent  operator  equation  actinq  on  the  Banach  space  Cn. 


Lemma  3.1:  y(x)  is  a  solution  of  BV(S)  if  and  only  if  it  is  a  solution  of  the  operator 
equation 


where 


with 


L  = 


Ly  =  F 

■  I 

"Vl2 

-V21 

I-K2(A22-A22) 

‘*2*23 

,  and 

_-£K3A31 

‘*3*32 

I 

1 

1 

e~1Y1<*,0)g*S) 

+  V, 

F  = 

Y2(''°>q2S) 

+  K  f 

2  2 

e“1Y3C,1)g5S)  +  V3 


(i^w)  (x)  3  e_1  J*  Y^ (x,s)w( s)  ds  , 

( K^w)  (*)  =  Iq  Y2(x,s)w(s)  ds  ,  and 

(x3w)  (x)  3  e-1  /*  Yj (x,s)w( s)  ds  . 


Proof:  Consider  the  initial-value  problem 

~^-w(x)  -  C(x)w(x)  =  h(x)  and  w(0)  =  h  , 
whose  fundamental  solution  matrix  Y(x,s)  satisfies  the  matrix  initial-value  problem 

3 

r —  Y(x,s)  -  C(x)Y( x , s)  -  0  and  Y(s,s)  *  1 
Sx 

The  variation  of  parameters  formula  states  that  the  solution  of  this  initial-value  problem 
admits  the  representation 
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w(x)  -  Y(x,0)n  +  /*  Y<x,s)h(s)  ds 


Suppose  then  that  y(x)  is  a  solution  of  BV(S).  Applying  the  variation  of  parameters 
formula  to  each  component  of  the  differential  equation  Ly  =*  f,  and  the  boundary- 
conditions  specified  by  B^s^y  «  g*s^,  shows  that  y(x)  is  also  a  solution  of  the 
operator  equation  Ly  -  F.  Conversely,  suppose  that  y(x)  is  a  solution  of  the  operator 
equation  Ly  =  F.  Then  y(x)  clearly  satisfies  the  boundary  conditions  8 *s'y  -  g*s'  and 
differentiation  of  each  component  of  the  operator  equation  shows  that  y(x)  is  also  a 
solution  of  Ly  «  f. 

D 

To  solve  the  integral  equation  Ly  *>  F  note  that  it  is  possible  to  compute,  for  all 
sufficiently  small  positive  e,  an  approximate  inverse  of  L. 


Lemma  3.2:  For  all  sufficiently  small  positive  e 


J 0 L  =*  I  -  CM 


where 
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l  3  32J 

[K2A21  1  K2R231  ' 


K1A12K2A23K3A31 


K2A23K3A31 


*3<A32K2A23VI)A31 


K1A12N2  K1(A12K2A21K1+I)A13 


*3*32*2 


K2A21K1A13 


*3*32*2*21*1*13 


e  k2[a22  +  A21K,A12  +  A,,K,A 
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'23  3ft32 


*225 


are  operators  on  Cn  which  are  bounded  independently  of  e . 


and 


Proof:  It  is  a  simple  matter  to  multiply  JQ  by  L  and  verify  the  product  is  indeed 
I  -  eM.  To  show  that  Jq  and  K  are  bounded  independently  of  e,  for  all  sufficiently 
small  positive  e,  one  uses  Theorem  2.3.  For  example,  to  prove  that  N2  is  bounded  first 
change  the  order  of  integration  in  the  double  integrals  yielding  (N2w)(x),  substitute  for 
the  fundamental  solution  matrices  Y^  and  Y^  the  expressions 
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Y.(x,s)  =  e  A  ]  ( x)  Y  (x,s) ,  and  Y  ( x ,s)  =  e  a”}(x)  |—  Y  (x,s)  , 

1  llaxi  3  33  d  X  3 

then  perform  an  integration  by  parts  and  recall  the  definition  of  A^t*)  before  applying 

Theorem  2.1. 

a 

Using  Theorem  2.4  we  can  now  obtain,  for  all  sufficiently  small  positive  e,  a 
convergent  generalized  power  series  expansion  for  L-\  the  fundamental  solution  matrix 
Y*s',  and  the  solution. 


Theorem  3.4:  For  all  sufficiently  small  positive  e,  L' 
generalized  power  series  expansion: 


L-1  =  [  I 
1=0 


exists  and  admits  the  following 


/  e  \ 

Corollary  3.5;  For  all  sufficiently  small  positive  e,  Y'  exists  and  admits  the 
following  generalized  power  series  expansion: 


Y<S>  -  [  l 


£=0 


,0) 


y2(*  ,0) 


e  1y;j(. 


,D 


Corollary  3.6:  For  all  sufficiently  small  positive  e,  BV(S)  has  a  unique  solution  y 
which  admits  the  following  generalized  power  series  expansion: 

Qo 

y  =  [  T  <em/1j0 
£«0 

As  shown  by  Theorem  2.4,  when  the  above  infinite  sums  are  truncated  after  the  term 
f.=N  the  resulting  approximations  to  L-1,  Y^S',  and  y  are  Q.(eN)  accurate  as  e  ♦  0  +  . 

When  the  solution  y  of  Ly  *  f  satisfies  Requirement  1.1  the  following  Corollary 
shows  that  y2<x)  remains  bounded  on  [0,1]  as  e  0+.  A  more  detailed  calculation, 


-1  (R) 

e  Y  ^  ( •  ,  0 )  q  ^  +  K^f  1 


Y2C,0,q'S)  ♦  V2 


e_1Y_  (•  .0)afS)  +  K  f  ' 
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such  as  that  provided  by  the  method  of  matched  asymptotic  expansions  [3],  is  needed  to  show 
that  when  the  solution  y  of  iy  »  f  satisfies  Requirement  1.1  both  y,(x)  and  V3(:  . 
remain  bounded  on  each  fixed  open  subinterval  of  [0,1]  as  e  +  0+. 

Corollary  3.7;  For  all  sufficiently  small  positive  e,  BV(S)  has  a  unique  solution  y 
which  admits  the  following  a  priori  bounds: 

"V  i  K  1  I*,*0''  +  'y2<o)l  +  e|y3n>!  +  « f ,»  +  * f 2* ,  +  « «3» >  » 

*  Y2 II  <  K  {e ly, (0) !  +  |y2<0)l  +  e  !  v3  { 1 )  I  +  If,!  +  Ifjl,  +  lf3D  .  and 

* Y3«  <  K  {e |y, (0 ) I  +  I Vj ( 0 ) I  +  lv3(1)l  +  lf,l  +  lf2l,  +  lf3D  • 

Here  the  norm  Ifl  of  a  function  f(x)  is  defined  as 

» f • 1  i  /glf(x)l  dx  . 
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4.  Solution  of  the  General  Problem 

By  combining  Theorem  2.2  and  Corollary  3.5  one  arrives  at  a  explicit  method  for 
determining  the  unigue  solvability  of  BV(G).  Recall  that  by  Theorem  2.2,  BV(G)  is 
uniquely  solvable  if  and  only  if  8<G)Y(S)  is  nonsingular.  To  arrive  at  reasonably 
general  conditions  under  which  the  solution  of  BV(G)  exists  and  satisfies  Requirement  1. 1 
let  us  follow  Harris  [6]  and  introduce  the  following  definition. 

Definition  4.1:  BV(G)  is  said  to  be  regular  if  and  only  if  it  admits  Assumptions  2.1  and 
the  matrix 

R-  lim+  8(GVS) 

e*0 

exists  and  is  nonsingular. 

Note  that  the  smoothness  requirement  of  Assumptions  2.1  allow  us  to  compute  R  using 
the  equivalent  formula 

" 

•"’V'O) 

R=  lim+  B<G)(I+€M)J0  Y2(*,0) 

£♦0 

e'V(*  ,1) 

L  J 

Theorem  4.2:  If  BV(G)  is  regular  then  for  all  sufficiently  small  positive  .  e  it  is 
uniquely  solvable  and  its  solution  y  satisfies  Requirement  1.1. 

Proof:  By  a  continuity  argument  the  nonsingularity  of  R  and  the  assumed  smoothness 

properties  of  0V(G)  guarantee  that  g(G)y<s>  ig  nonsingular  for  all  sufficiently  small 

positive  c.  Therefore,  by  Theorem  2.2,  BV(G)  is  uniquely  solvable  for  all  sufficiently 

small  positive  E.  To  verify  that  the  solution  of  BV(G)  satisfies  Requirement  1.1  let  us 

observe  that  it  admits  the  representation 

e'V^O,.  Key^O)}  + 

YJO,*  ){  y2(0)}  +  K2f2 

E_,Y3(*,1){Ey3(1)}  +  K3f3 


y  -  (  l  <eM)*]J0 

t-0 
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The  smoothness  requirements  In  Assumptions  2.1  and  the  regularity  of  BV(G)  assure  us  that 
the  limit 


lim+ 

e+0 


ey^O) 

(  /  00 

Vi 

y2(0) 

-li»+  (b<GVS>  -1  g(G)-B<G)[I  (£M)^]  J. 

£■*•0  t- 0 

K„f  - 
2  2 

ey3(1) 

^3  j 

) 

exists  and  is  finite. 
□ 
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5.  Conclusions 


Similar  techlques  can  be  used  to  establish  the  stability  and  consistency  properties  of 
finite  difference  approximations  to  BV(v)  for  v  -  G,  S  based  on  Euler's  method.  This 
finite  differenca  scheme  [3,10],  with  step  size  h,  uses  the  backward  Euler  scheme  to 
approximate  the  differential  equation  for  y^,  the  box  scheme  of  Keller  (centered  Euler) 
to  approximate  the  differential  equation  for  y2,  and  the  forward  Euler  scheme  to 
approximate  the  differential  equation  for  y^.  This  scheme  is  Q.(h^)  accurate  as  h  ♦  0+, 
uniformly  in  e,  in  each  fixed  open  subinterval  of  [0,1]. 

The  a  priori  estimates  of  Corollary  3.7  can  also  be  used,  as  in  [3] ,  to  show  that  the 
method  of  matched  asymptotic  expansions  does  indeed  yield  an  asymptotically  correct 
estimate  of  the  solution  of  BV(G)  as  e  +  0+.  These  a  priori  estimates  can  also  be  used 
to  prove  the  existence  of  solutions  of  a  class  of  nonlinear  singularly  perturbed  boundary- 
value  problems  when  the  method  of  matched  asymptotic  expansions  is  coupled  with  the  Newton- 
Kantorovich  theorem. 

More  general  boundary-value  problems  of  the  form 

')d_jV*>]  fcn<*>  Ci2<xf|r*i(xn  ffl(x>i 
ei]  d*  [*2<*)J  '  fc21(x>  C22<x)J(z2<x)J  "  J2<x)J  '  Z  “  9 

can  also  be  analyzed  by  this  technique  under  the  restriction  that  the  eigenvalues  of 
C22(x)  have  real  parts  whose  absolute  value  is  bounded  away  from  zero.  For  such  problems 
a  preliminary  block  diagonalization  of  C22(x)  must  be  performed  before  the  results  of 
sections  3  and  4  can  be  applied.  This  block  diagonalization  of  C22(x)  Involves  the 
existence  of  a  smooth  nonsingular  matrix  U(x)  with  the  property  that 

C22(X\ 

where  the  real  parts  of  the  eigenvalues  of  Cjjtx)  ar®  negative  (positive).  The 

existence  of  such  a  similarity  transformation  is  established  in  [3,5]. 

These  results  summarize  part  of  a  thesis  [3]  written  under  the  direction  of  Dr.  H.  B. 
Keller.  I  would  like  to  thank  Dr.  Keller  for  suggesting  this  problem  to  me  and  for  his 
supervision  of  my  studies  at  the  California  Institute  of  Technology. 
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